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Orbital  resonances  are  widespread  throughout  the  solar  system,  Itiey 
appear  mainly  between  the  satellites  of  the  Jovian  planets,  but  also 
are  found  among  the  planets  themselves  as  In  the  case  of  the  3/2  res- 
onance between  Pluto  and  Neptune,  This  project  undertook  to  examine 
the  conditions  out  of  which  a resonance  could  be  established. 

Of  particular  interest  was  the  3/2  resonance  of  Pluto  and  Neptune. 
This  resonance  has  been  firmly  established  by  a number  of  researchers. 

This  project  also  relates  to  the  hypothesis  that  Pluto  might  have  origin- 
ated as  a moon  of  Neptune  that  esoaped  and  entered  an  orbit  independent 
of  Neptune,  The  last  major  objection  to  this  is  that  Neptune  and  Pluto 
never  dome  close  to  one  another  since  they  are  in  resonance.  The  lack 
of  a close  encounter  between  Pluto  and  Neptune  is  due  to  the  resonance 
and  not  a consequence  of  it.  Through  the  use  of  numerical  integration 
this  paper  establishes  that  it  la  possible  that  if  Pluto  escaped  from 
Neptune  then  it  could  go  into  an  orbit  and  a resonance  very  similar 
to  that  which  exists  today.  Once  that  resonance  is  established  it  is 
impossible  to  regain  the  initial  conditions  due  to  the  fact  that  a 
resonance  is  self-perpetuating.  The  integrations  were  carried  out  for 
a time  equal  to  over  one  million  years  in  a model  solar  system  containing 
all  the  planets  from  Jupiter  outwards.  This  paper  removes  that  last 
major  objection  to  the  hypothesis  that  Pluto  may  have  originated  as  a 
moon  of  Neptune, 

An  additional  result  was  the  development  of  a system  for  identifying 
resonances  through  the  use  of  a rotating  coordinate  system. 
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Introduction  and  Proposal 

Throughout  the  solar  system  there  are  many  examples  of  orbital 
resonance.  An  orbital  resonance  is  defined  as  any  system  in  which  two 
or  More  satellites  are  orbiting  a primary  such  that  the  ratio  of  thtlr 
periods  is  a ratio  of  smII  whole  numbers.  Periods  of  this  type  ax* 
said  to  be  commensurate  and  the  resonance  is  referred  to  as  a commensur- 
abillty.  When  discussing  satellites  and  primaries  included  also  are  the 
motions  of  the  planets  around  the  Sun. 

One  of  the  best  examples  of  a simple  two  satellite  resonance  exists 
between  two  moons  of  Saturn,  Titan  and  Hyperion.  Here  a large  massive 
moon.  Titan,  moves  in  a nearly  circular  orbit.  Hyperion,  however,  is 
smaller  and  moves  in  an  orbit  of  large  eccentricity.  IMs  system  has 
been  studied  extensively  and  it  has  been  found  that  the  ratio  of  the 
periods  is  equal  to  4/3,  That  is,  if  n^  and  are  the  periods  of  Titan 
and  Hyperion,  respectively,  then  the  relation  4n^  - 3n^  * 0 is  satisfied 
to  a very  high  degree.  It  has  also  been  determined  that  the  point  of 
conjunction  or  the  point  of  closest  approach  between  the  moons  oscillates 
ahout  the  apocenter  of  the  orbit  of  Kyp*rlon  with  a maximum  amplitude  of 
36°  and  a period  of  approx iamately  2 years,1  In  ad lltion  to  these 
oscillations  the  line  of  apsides  about  which  the  oscillations  occur 
processes  in  a retrograde  sense,  opooslte  to  the  motion  of  the  bodies, 
at  a rate  of  about  t9°  per  year.  This  is  also  a result  of  the  resonance 
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The  Titan-Hyperion  resonance  is  a classic  example  of  a large  category 
of  orbital  resonances  termed  large  eccentricity  resonances.  The  mechanism 
for  this  resonance  is  relatively  simple.  Suppose  that  a primary  has  two 
satellites  in  orbit  around  it.  If  the  moons,  m^  and  m^,  are  such  that 
*^»  «2  and  they  are  in  coplanar  orbits  such  that  is  in  an  inner 
circular  orbit  and  m2  is  in  an  outer  eccentric  orbit  then  the  situation 
shown  in  Figure  1 exists.  If  initially  conjunction  occurs  at  point  A 
then  shortly  before  conjunction  the  two  moons  are  in  positions  labeled 
#1.  The  gravitational  force  between  them  can  be  resolved  into  two  com- 
ponents. one  normal  inwards  towards  the  primary  and  one  tangential, 
opposite  to  the  motion  of  n^.  Similar  forces  act  on  but  are  negligible 
due  to  the  large  mass  of  in  comparison  with  m2.  The  net  effect  of 
these  forces  is  to  remove  angular  momentum  from  the  small  moon.  Once 
past  conjunction  there  exists  a time  when  the  satellites  are  in  position 
#2  such  that  ©^  • ©2>  Here  also  the  force  has  a normal  and  a tangential 
component.  After  conjunction,  however,  the  tangential  component  adds 
angular  momentum  to  the  smaller  satellite.  Since,  at  point  A the  orbits 
are  diverging,  the  distance  between  and  ir^  is  greater  at  point  *2 
than  it  was  at  point  #1,  Therefore  the  force  is  less  and  the  tangential 
component  does  not  add  back  as  much  angular  momentum  as  was  lost  before 
conjunction.  The  net  overall  effect  after  one  cycle  is  that  m2  loses 
angular  momentum.  To  stay  in  orbit  it  must  speed  up.  If  the  orbits 
are  commensurate  then  the  next  time  conjunction  occurs  it  will  be  at 
some  point  slightly  advanced  along  the  orbit. 

The  situation  is  reversed  if  the  conjunction  occurs  at  B.  Now  the 
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additive  tangential  component  dominates  and  the  smaller  moon  gains 


angular  momentum  therefore  slowing  down.  The  successive  conjunctions 
will  be  displaced  back  towards  the  line  of  apsides.  If  conjunction 
occurs  on  the  line  of  apsides  then  the  net  effect  of  the  components  can- 
cel and  the  conjunction  point  is  stable*  We  therefore  see  that  the 
resonance  keeps  the  orbits  stable  and  that  the  point  of  conjunction  tends 
to  oscillate  about  the  apocenter.  The  other  end  of  the  line  of  apsides, 
the  perioenter.  is  also  an  equilibrium  point  but  It  is  an  unstable 
equilibrium.  Any  deviation  away  froe  this  point  causes  the  system  to 
begin  to  oscillate  about  the  apocenter.  The  oscillations  help  preserve 
the  comaensurabllity  of  the  system. 

This  is  exactly  the  situation  of  Titan  and  Hyparlon.  They  are  in 
oommensurmte  orbits  and  show  an  oscillation  of  the  conjunction  about  the 
apocenter  of  36°  over  a 2 year  period.  Also  present  in  the  Titan-Hyperlon 
system  is  the  retrograde  motion  of  the  line  of  apsides.  This  can  be 
explained  with  large  eccentricity  resonances  also. 

In  a system  in  which  the  point  of  conjunction  is  oscillating  it 
will  at  some  time  fall  on  the  apocenter.  The  force  is  now  directed 
inwards  the  primary  axolttletlfr.  ‘nils  force  tends  to  pull  the  outer 
moon  in  tasuti  the  primary  and  it  can  be  thought  of  as  having  traveled 
on  an  orbit  slightly  inside  the  unperturbed  orbit.  Extending  this  orbit 
around  we  find  that  will  reach  its  closest  point  of  approach  slightly 
sooner  than  it  would  have  otherwise.  It  appears  then  that  the  line  of 
apsides  has  rotated  in  a direction  opposite  to  that  of  the  satellite 
notion  or  thus  in  a retrograde  manner.  This  is  exactly  the  case  in  the 
Titan-Hyperlon  reaoaance.  The  retrograde  motion  completely  dominates 
any  other  pro-r^de  motions  and  the  line  of  apsides  rotates  approxiamately 
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19°  P*r  year. 

The  Tltan-Kyoerlon  case  is  only  one  example  of  many  systems  that 
exhibit  these  resonance  phenomenon.  Resonances  have  been  detected  between 
most  of  the  moons  of  the  outer  planet,  in  the  rings  of  Saturn,  within 
the  asteroid  belt,  and  even  among  the  planets  themselves. 

Of  primary  interest  to  this  report  are  the  orbits  of  the  planets 
Pluto  and  Neptune,  There  are  a great  many  facets  about  the  orbits 
which  point  out  that  some  kind  of  resonance  exists  between  the  two. 

The  first  factor  is  that  the  orbital  periods  are  approxlamately  commen- 
surate in  a ratio  of  3 to  2,  This  leads  us  to  suspect  a 3/2  resonance. 

In  addition,  this  resonance  is  of  the  large  eccentricity  type.  The 
system  has  all  the  correct  initial  conditional  the  mass  of  Neptune  is 
on  the  order  of  100  times  greater  than  that  of  Pluto,  satisfying  the 
m^»a2  oondltlon.  Pluto's  orbital  *ocentrlcity  is  the  largest  of  the 
planets)  so  large  th*t  Pluto's  orbit  crosses  inside  of  Neptune's  orbit. 
Neptune,  on  the  other  hand,  has  one  of  the  smallest  eccentricities  putting 
it  in  essentially  a circular  orbit.  The  system  has  all  the  necessary 
conditions  for  a large  eccentricity  resonance. 

If  we  look  at  the  Pluto-Neptune  system  today  we  see  the  same  effects 
that  are  present  in  the  Titan-Hyperion  system  that  characterize  an  orbital 
resonance.  Conjunctions  oetween  the  two  planets  oscillate  about  the 
apocenter  with  the  apocenter  as  the  stable  equilibrium  point.  Also  the 
line  of  apsides  of  Pluto  processes  in  a retrograde  motion.  All  of  these 
similarities  with  the  Titan*Hyperion  system  have  led  us  to  believe  that 
today  Pluto  and  Neptune  exist  in  a 3/2  large  eccentricity  resonance. 

In  addition  to  the  features  mentioned  above,  Pluto  exhibits  a third 
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phenomenon  possibly  caused  by  the  eommensurability . This  is  that  the 
resonance  appears  to  keep  the  conjunctions  near  90°  from  the  mutual 
orbital  nodes,  this  eliminates  the  possibility  of  a close  enoounter 
or  oollision  between  Neptune  and  Pluto. 

The  major  question  that  arises  1st  how  did  this  resonanoe  wet  started? 
One  theory  is  that  Pluto  is  an  escaped  moon  of  Neptune.  There  are  many 
faotors  that  seem  to  make  this  at  least  possible.  The  first  is  that 
the  orbit  of  Pluto  comes  dose  to  and  actually  inside  that  of  Neptune. 

A second  fact  is  that  the  mass  of  Pluto  is  low  enough  for  it  to  be  a 
moon  of  Neptune.  In  fact*  the  mass  of  Pluto  is  of  the  same  order  of 
magnitude  as  one  of  Neptune* s current  moons.  Triton. 

If  we  theorise  that  Pluto  was  at  one  time  in  orbit  around  Neptune 
then,  due  to  the  tidal  effects  between  the  two.  we  expect  Pluto's  rot* 
atlon  to  have  slowed  to  the  point  where  it  would  have  a period  of 
rotation  equal  to  its  period  of  revolution  about  Neptune.  In  other 
words,  Pluto  would  always  keep. the  same  side  towards  Neptune  much  as 
our  swn  Noon  does  with  resoect  to  the  Sarth.  After  ejection  from  orbit 
Pluto  might  keep  the  same  oeriod  of  rotation.  Indeed  measured  values 
for  the  period  are  very  close  to  what  we  would  expect  for  a large  satellite 
of  Neptune,  approxlamately  5,5  x 10^  seconds.  If  we  use  Kepler's  Third 

i 

Law  which  relates  the  period  to  the  distance  from  the  primary  to  compute 

g 

a value  for  an  orbital  radius  around  Neptune  we  obtain  3.$  x 10  meters. 
This  is  very  close  to  the  orbital  radius  of  Neptune's  moon, 

A third  piece  of  evidence  that  Pluto  might  have  evolved  out  of  the 
Neptune  system  of  moons  is  that  it  appears  that  the  inclination  of  Pluto's 
orbit  to  the  "eaHptie  (approxlamately  17°)  is  about  the  same  inclination 
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as  the  orbits  of  Neptune's  moons,  Pluto's  large  inclination,  as  con- 
trasted with  that  of  the  other  planets  wnlch  have  Inclinations  of  at 
nost  3*5°  (Venus  - 3-39°),  makes  this  seem  more  than  a coincidence. 

The  last  and  possibly  most  conclusive  piece  of  evidence  In  favor 
of  Pluto  being  an  escaped  moon  of  Neptune  is  the  Irregular  structure  of 
the  orbits  of  Neptune's  present  moons,  Neptune  nas  two  moonsi  Triton, 
which  Is  slightly  smaller  than  Pluto,  and  Nereid,  which  resembles  a 
stray  asteroid.  The  Irregularity  is  that  Triton,  although  in  a nearly 
circular  orbit,  revolves  In  a retrograde  manner.  This  In  itself  is  not 
too  unusual  but  then  we  notice  that  Triton  Is  the  only  massive  sa'^Jlite 
In  the  solar  system  to  exhibit  this  retrograde  motioni  all  the  ether 
satellites  with  retrograde  motion  are  small  enough  to  be  captured  asteroids, 

I 

If  we  assume  that  Triton  and  Neptune  were  created  at  the  same  time, 

which  is  not  unreasonable,  then  something  must  have  happened  to  cause 

Trltdn  to  turn  around  In  Its  orbit.  That  something  is  thought  to  be 

2 

a close  or  near  collision  with  Pluto,  If  at  one  time  Pluto  and  Triton 
were  both  In  orbit  around  Neptune  then  the  situation  exists  whereby  the 
two  could  have  an  interaction  which  boosted  Pluto's  speed  so  It  could 
escape  from  Neptune  and  go  into  orbit  around  the  Sun  on  its  own.  A 
corollary  to  this  is  the  fact  that  Pluto  resembles  a noon,  small  and 
icy,  more  than  It  resembles  the  gas  giants  of  Jupiter,  Saturn,  Uranus, 
and  Neptune  which  populate  the  outer  solar  system. 

However,  one  difficulty  results  In  this  theory.  If  Pluto  did 
originate  as  a moon  of  Neptune  we  would  expect  that  Pluto  and  Neptune 
would  ocasi&aally  come  close  to  each  other  due  to  the  cyclical  nature 
of  the  orbits.  They  do  not!  In  fact  the  orbits  of  Pluto  and  Neptune 
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never  do  actually  Intersect  and  the  planets  never  get  close  to  one  <■ 
another,  an  observation  which  has  been  used  against  the  Neptune-moon 
hypothesis  for  the  origin  of  Pluto.  The  reason  behind  this  may  be  due 
to  the  establishment  of  a resonance. 

The  resonance  between  Neptune  and  Pluto  has  been  established 
through  numerical  integrations  and  orbital  observations  to  be  stable 
In  a 3/2  ccifimensurability,^  Therefore,  If  the  present  situation  is 
taken  and  Integrated  both  backward  and  forward  in  time  the  orbits  will 
not  change  significantly  because  the  system  Is  locked  into  the  resonance. 
The  problem  is  that,  under  resonant  conditions.  Integrating  backwards 
will  never  bring  us  to  the  initial  conditions  of  the  escape  that  set 
up  the  resonance.  The  analogy  may  be  made  to  the  case  of  an  object 
falling  through  the  atmosphere  with  terminal  velocity  1 observation  of 
the  final  state  provides  no  information  about  the  intial  conditions, 
fewewarif  we  start  with  the  moment  of  escape  and  integrate  forward  in 
time,  then  it  might  be  determined  what  conditions  are  necessary  for  the 
establishment  of  the  present  3/2  resonance.  Such  is  the  intent  of  this 
paper. 
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In  order  to  Integrate,  over  time,  the  position  end  velocity  of 


spatial  bodies,  we  muet  first  have  a coordinate  system  for  a reference 


The  elapleat  and  easiest  to  work  with  Is  a Cartesian  reference  frane 
centered  on  the  Sun.  This  le  illustrated  In  Figure  2,  The  point  (x,y,a) 


represents  the  position  of  Pluto  at  soma  time  t.  The  acceleration  due 


to  the  Sun  la  written  aai 


G « Universal 

C.ravi  tationnl 
constant 


In  this  fora  the  equations  are  relatively  easy  to  work  with.  Per 


turbatlona  due  to  Neptune  and  the  other  planets  are  added  with  the  add 


ltlon  of  similar  terms  of  the  same  form  except  that  the  distances  of 


the  perturbing  planet  relative  to  Pluto  are  substituted  for  x,  y,  and 
a and  the  perturbing  planets  mass  is  substituted  for  the  Sun's  mass. 
However,  the  Interpretation  of  the  results  la  difficult  If  we  are  looking 


for  changes  In  one  orbit  relative  to  another.  It  would  be  beneficial 


if  we  oould  fix  one  rlanet  and  observe  only  Pluto's  motion  relative  to 


This  can  be  done  using  a rotating  coordinate  frame  fixed  with 


mm 


‘V1*1'  r*WWHti 


acos(wt)  - bsin(wt)  -2w(asln(wt)+bcos(wt))-w2(acos(wt)-bsln(wt) ) 

“ -^CM  ^ yz  (acos(wt)  - bsin(wt)) 


iln(wt)  ♦ bcos(wt)  +2w(acos(wt)-bsln(wt))  -w2(asln(wt)-*-bcos(wt)) 


(aff+b^+z2)  V2  (asln(wt)  ♦ bsln(wt)), 


These  are  Eqs.  (6a)  and  (6b), 


To  simplify  these  we  will  leave  out  the  z term  as  the  z coordinate  trans- 
forms unchanged.  We  now  use  a mathematical  trick  of  multiplying  the  first 
by  cos(wt)  and  the  second  by  sln(wt)  and  then  adding  the  two.  Letting 


s • sln(wt)  and  c - cos(wt)i 

ac2  -bsc  -2w(asc  + be2)  - w^ac"  - bsc) 


GK 

r^Ti 


(ac2  - bsc) 


as2  ♦bsc  +2w(asc  - bs2)  - w2(as2  + bsc)  - - ^ (as2  ♦ bee) 
Adding  and  using  s2  + c2  - sln2(wt)  + cos2(wt)  - 1 


a - 2wb  - w2a  - -(JT^Tp/2  <a)  • (?a) 

Performing  a similar  operation  againi  only  multiplying  the  first  by 
sln(wt)  and  the  second  by  cos(wt)  and  then  subtracting  the  two.  After 
exoandlng  and  collecting  all  terms  we  are  left  withi 

b + 2wa  - w2b  - (b).  (?b) 

Equations  (7a)  and  (7b)  can  be  written  to  give  the  accelerations 
in  the  rotating  (A, 3)  coordinate  frame i 

» " " TIzfep/2  (a)  * ^a  ♦ 2wl>  (8a) 

b - -(;2Vb?p/2  (b)  * w2b  - 2wa  . (8b) 
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The  equations  are  what  Is  expected  In  a rotating  frame  with  the 
first  term  the  gravitational  force,  the  second,  the  centrifugal  force, 
and  the  third,  the  corlolls  term. 


Now  that  the  accelerations  In  (A, 3)  have  been  determined  the  position 
and  the  velocity  transformations  must  be  worked  out.  The  position 
transformations  have  already  been  given  In  Sq,  (2) i 

x ■ aco8(wt)  - bsln(wt) 
y - asln(wt)  +■  bcos(wt)  , 

Writing  this  as  a matrix  equation  we  havei 


[;]•  [: 


cos(vt)  -sin(wt) 
sin(wt)  cos( wt ) 


Defining  T as  the  transformation  matrix  such  that 


[’cos(wt)  -sin(wt) 

sln(wt)  cos(wt) 


Then  to  go  from  the  non-rotating  (X,Y)  to  the  rotating  (A, 3),  we  need 


only  to  Invert  the  matrix  and  apoly  it  to  (a,b) 

cos(wt)  8ln(wt) 

T ■ 

— -sln(wt)  eos(wt) 


Thus, 


a -EL1  x 

b y 


’cos(wt) 

-sin(wt) 


sln(wt) 

cos(wt) 


a • xcos(wt)  + ysin(wt) 
b - -xsln(wt)  + ycos(wt)  . (9) 

To  obtain  the  velocity  transformations  we  go  back  to  Sqs.  (3a)  and 

(3b)  which  result  in,  after  rearrange  Inn 
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x - acos(wt)  - bsln(wt)  - wy 
y - asln(wt)  bcos(wt)  + wx  . (10) 

Eqs,  (3)  can  also  be  written  asi 

x - (a  - wb)cos(wt)  - (b  + wa)sin(wt) 
y ■ (a  - wb)sin(wt)  + (b  + wa)cos(wt)  , (11) 

These  two  sets  of  equations  give  us  the  transformation  matrix  between 

velocities  in  the  two  coordinate  systems i 


x m cos(wt)  -sin(wt)  a-wb 
y sin(wt)  cos(wt)  b*wa 

U ^ 

and  thus, 


♦ ■ T { 

y — b+-wa 


Rewritting  £q,  (10),  we  obtain 


x + wy 
y - wx 


Inverting  this  gives i 


co8(wt)  -sin(wt) 
sin(wt)  cos(wt) 


a |cos(wt)  sin(wt)  x + wy' 

b -sin(wt)  cos(wt)  y - wx 


i • r.1 


x + wy 
y - wx 


Now  we  have  a complete  set  of  equations  and  coordinate  transformations 
between  fixed  (X.Y-)  coordinate  frame  and  rotating  (A,R)  frame. 


Summarising! 


Non-rotatlr 


..  . GMx 

X m(*Gy**i 


Rotatlr 


•%  GMa  2 • 

a " ~(^Vb^^)3/2  a 


b - -/ 


+w2a  -2wi 


Rotating 
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a*+b*+z< 


To  transform  from  non-rotating  to  rotalng  coordinates  knowing 

x»y**»x,y,ii 

a - xcos(wt)  + ysin(wt) 
b - -X8in(wt)  +ycos(wt) 

Z • Z 

a - (x  + wy)cos(wt)  + (y  - wx)sin(wt) 

b - -(x  + wy)sin(wt)  + (y  - wx)cos(wt) 

• « 
s * z 

Similarly  to  go  from  rotating  to  non-rotating  reference  frames  knowing 

. * • * 
ftfbf  E|Af  Df  Z I 

x - acos(wt)  - bsln(wt) 
y - asin(wt)  + bcos(wt) 
z - z 


x - (a  - wb)cos(wt)  - (b  + wa)sin(wt) 
y - (a  - wb)sln(wt)  + (Ij  +■  wa)cos(wt) 


9 9 

S • Z 


In  the  extreme  case  where  wt-O,  which  will  appear  later  when  we 
wish  to  start  on  the  x-axis,  the  transformations  from  the  fixed  to 
rotating  reduce  to  1 


a - x 
b - y 


z - z 


a - x + wy 
b ■ y - wx 
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Recognition  of  Resonances 

By  transforming  to  the  rotating  coordinate  system  three  advantages 
are  gained.  Firstly,  the  motion  of  one  of  the  planets  Is  eliminated 
from  the  picture  so  that  now  only  the  motion  of  the  perturbed  planet 
relative  to  the  perturbing  planet  is  observed.  This  is  exactly  what  is 
needed  for  a study  of  resonant  motion.  Secondly,  by  fixing  one  planet 
the  calculations  become  more  accurate  over  long  time  periods  because  any 
numerical  integration  of  perturbations  is  now  referenced  to  a fixed 
non-calculated  point  rather  than  to  a moving  calculated  point.  And 
thirdly,  by  going  to  a rotaing  coordinate  system  the  graphical  display  of 
the  perturbed  planet's  motion  provides  an  easy  method  for  determining 
the  exact  resonance  or  commens urabi 11 ty  that  the  planets  are  in.  Thus, 
one  of  the  first  results  obtained  was  the  development  of  this  technique 
for  identying  resonant  orbits. 

In  order  to  identify  a resonant  orbit  it  is  necessary  merely  to 
plot  the  motion  of  one  planet  relative  to  another  in  a rotating  coordinate 
frame.  Such  a plot  is  shown  in  Figure  4.  Remembering  that  the  interior 
planet  is  always  fixed  6n  the  positive  x-axis  it  can  be  shown  that  the 
loops  are  the  oerihelion  passage  of  the  outer  planet.  The  number  of 
loops  is  characteristic  of  the  particular  orbital  resonance.  If  the 
resonance  type  (i.e.  the  ratio  of  the  periods)  is  expressed  as  a fraction, 
for  example  3/2  or  4/3,  the  number  of  loops  per  complete  pattern  is  the 
denominator  of  the  fraction  or  ratio.  The  relative  orbital  period  df  the 
perturbing  planet,  indicated  by  the  numerator  of  the  fraction  can  be 
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determined  by  observing  bow  many  orbits  It  takes  for  the  pattern  to 
oomplete  Itself.  Thus,  if  It  takes  3 orbits  to  complete  a two  loop 
pattern,  a 3/2  resonance  is  present.  Ibis  Is  the  case  In  Figure  4, 

An  easier,  although  more  tiae  consuming  method  of  determining  a pattern 
is  to  plot  the  x coordinate  of  the  motion  versus  time.  The  time 
necessary  to  complete  one  cycle  is  the  same  as  the  time  needed  to  com- 
plete pattern. 

By  way  of  illustration  the  following  are  a number  of  different 
resonant  patterns  with  their  x vs.  t plots. 


>1 
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FIGURE  8B,  X vs,  T plot  for  a 2/1  Resonance  with  cycle  complete  at  T 


All  the  patterns  shown  are  in  their  stable  non-oscillating  config- 
urations. In  an  actual  oscillating  resonant  display  the  only  difference 
is  that  the  looped  pattern  itself  oscillates  a rout  a point.  The  ampli- 
tude of  these  oscillations  is  directly  related  to  the  amplitude  of  the 
actual  conjunction  oscillations.  Figure  10  shows  this  for  an  oscillating 
3/2  resonance.  Note  that  the  point  of  stable  equi librium  lies  on  the 
positive  y-axis.  This  corresponds  to  the  aphelion  of  the  outer  eccentric 
orbit. 

A third  feature  of  a stable  resonant  pattern  is  seen  when  a three 
dimensional  picture  of  the  resonance  pattern  between  inclined  orbits 
such  as  those  that  exist  between  Pluto  and  Neptune  is  produced.  Since 
the  perihelion  lies  in  the  positive  z region  of  space  the  loops  also 
will  be  in  the  positive  z region  of  space.  This  is  apparent  in  Figure  11. 
Conversely  when  the  perihelion  lies  at  the  node  the  loops  degenerate 
into  curves  that  only  appear  as  loops  in  the  x-y  projection.  The  three 
dimensional  pictures  shown  here  were  generated  using  the  Evans-Sutherland 
Picture  System  in  conjunction  with  the  PDP-11  computer  and  plotted  on 
a XYnetics  flatbed  plotter.  The  two  dimensional  plots  were  obtained 
using  the  Tektroni*  4051  "smart”  graphics  terminals. 


FIGURE  10,  Oscillating  3/2  Resonance  In  Rotating  Coordinate  System 
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The  present  configuration  of  the  Pluto-Neptune  resonance  is  a 3/2 

commensurability  oscillating  about  a stable  equilibrium  point  with  a 

period  of  20,000  years  or  about  115  orbits  of  Neptune,  In  addition, 

there  are  a number  of  slower  oscillations,  the  most  noteworthy  of  jthich 

is  a 4 million  year  or  24,000  orbit  oscillation  of  the  argument  of  the 
U 

perihelion.  The  goal  then  is  a formation  process  starting  with  Pluto 
escaping  from  Neptune  and  eventually  reaching  something  close  to  the 
present  configuration  cited  above.  This  is  done  through  the  use  of 
several  computer  programs  to  numerically  integrate  the  positions  of 
the  planets  of  long  time  periods. 

The  computer  programs  used  to  perform  the  integrations  were  all 
written  in  BASIC  and  performed  on  the  U.  S.  Naval  Academy's  Dartmouth 
Time-Sharing  System  (DTSS),  Listings  of  all  programs  vised  are  given 
in  Appendices  II-  tf.  The  programs  were  designed  for  versatility  and 
ease  of  adaption  to  new  factors.  The  programs  consist  of  four  main 
parts.  First,  an  Input  section  where  variables  initialized.  Second, 
the  gravitational  force  of  the  Sun  and  any  perturbing  planets  are  laid 
out  in  the  x,  y,  and  z directions.  The  third  and  main  body  of  the  pro* 
grams  consists  of  a fourth  order  Rv&e-Kutta  numerical  integration  scheme 
where  the  forces  derived  in  the  second  section  are  evaluated  at  each 
of  the  four  steps  and  from  this  new  positions  and  velocities  of  Pluto 
are  computed.  The  fourth  and  last  section  is  made  up  of  various  output 
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statements  for  the  tabulation  of  the  results. 

Using  this  type  of  program  it  became  extremely  easy  to  add  or  subtract 
additional  planets  from  the  general  model.  Each  planet  appears  as  an 
additional  perturbation  to  the  planet's  motion.  Proceeding  in  this 
manner  the  planets  could  be  added  one  by  one  and  the  Influence  each  had 
on  the  developing  resonance  could  be  studied.  In  the  end  the  entire 
outer  flair  system  was  modeled. 

lb  overcome  any  Droblem  in  dealing  with  the  large  distances  con- 
cerned the  entire  solar  system  was  scaled  so  that  one  unit  of  distance 
became  the  Sun-Neptune  di stance}  one  unit  of  mass  is  equal  to  a solar 
massi  and  one  unit  of  time  was  taken  as  Neptune's  period  of  revolution 
about  the  Sun.  In  this  scaled  version  of  the  solar  system,  Newton's 
gravitational  constant,  G,  is  equal  to  the  orbital  speed  o^  Neptune 

is  2tn  and  thus  the  positions  and  velocities  of  the  other  planets  must 
be  sealed  accordingly. 

The  one  problem  that  develops  when  we  scale  the  solar  system  to 
Neptune's  level  is  that  the  period  of  the  revolutioiraf  Jupiter  and  Saturn 
are  so  much  shorter  than  Neptune's  that  during  one  integration  step 
Jupiter  and  Saturn  move  a fair  distance  along  their  orbits.  On  the 
other  hand  if  we  made  the  step  sise  smaller  to  limit  the  distance  Jupiter 
moves  then  the  run  times  for  the  program  exceeds  practical  limits,  A 
trade-off  must  be  made.  It  seems  that  a step  sise  of  0.01  orbits  of 
Neptune  is  the  best.  At  this  rate  it  takes  12.8  seoonds  of  CPU  (central 
processor  time)  to  complete  one  orbit  of  neptune,  168  years*  in  the 
model  solar  system  with  all  the  planets  past  Mars  present. 
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Assuming  Pluto  was  once  in  orbit  around  Neptune,  it  is  advantageous 

to  know  some  of  the  particulars  of  that  orbit.  Through  photometric 

studies  it  is  possible  to  determine  the  rotational  period  of  pluto  with 

5 

great  accuracy.  If  we  assume  that  none  of  Pluto’s  rotational  energy 
was  lost  in  any  escape  then  this  period  should  give  an  approxlamate 
value  for  Pluto’s  period  of  revolution.  Because  the  mass  differential 
between  Neptune  and  Pluto  is  so  great  it  is  reasonable  to  assume  that 
it  rotated  once  in  every  revolution  about  Neptune  much  the  smae  as  our 
Moon  does  around  Earth,  always  keeping  the  same  side  towards  the  Earth. 

Knowing  this  period,  P,  Kepler’s  Third  Law  can  be  used  to  compute  an 
orbital  radius  from  Neptune,  a, 

.3  - 2i£-n)  # 

Q 

This  value  is  a - J.8  x 10  m.  The  present  a for  Neptune's  other  massive 

Q 

moon,  Triton,  is  3.5  x 10  m.  It  is  then  conceivable  that  a near  collision 
between  the  two  resulted  in  reversing  Triton's  motion  and  ejecting  Pluto 
from  orbit.  See  Figure  12,  The  research  carried  out  here  does  not 
concern  itself  with  the  actual  collision  that  produced  the  escape  but 
only  the  result  of  such  an  escape  and  whether  it  was  possible  for  such  an 
escape  to  evolve  into  a 3/2  resonance  with  Neptune. 

A wide  range  of  initial  escapes  were  tried.  The  results  were  obtained 
first  from  a Meptune  centered  system  in  which  rotation  around  the  Sun 
was  added,  (See  program  NEPTUNE  in  Appendix  II)  This  introduces  a 
small  Coriolis  term  in  additon  to  the  Sun's  Gravity  into  the  force 
equations.  At  some  point  far  enough  away  to  avoid  scaling  problems  the 
coordinates  and  velocities  of  Pluto  were  transformed  to  the  Sun  centered 
reference  frame  and  allowed  to  continue  t*  its  orbit  around  the  Sun, 
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( See  program  C0NV2  in  Appendix  IY  for  conversion  factors).  By  scaling 
problems  we  mean  that  if  Pluto’s  position  is  transformed  into  the  Sun 
system  too  soon  Its  position  will  be  lost  in  the  significant  figures  that 
the  computer  can  handle.  Therefore,  care  must  be  taken  to  allow  Pluto 
to  travel  far  enough  away  from  Neptune  so  that  there  is  a significant 
separation  between  the  two  in  our  scaled  solar  system. 

The  results  indicate  that  Pluto  will  indeed  escape  if  given  even  a 
sizeable  fraction  of  the  necessary  escape  speed.  A resonance  between 
Pluto  and  Neptune  will  be  formed  in  a majority  of  the  cases.  These  range 
from  an  unstable  8/?  to  a very  stable  4/3  resonance.  See  Appendix  I for 
a discussion  of  relative  orbital  stability. 

The  3/2  resonance  will  be  established  if  a close  interaction  between 
Pluto  and  Triton  leaves  Pluto  with  near  escape  speed  heading  in  a direction 
that  is  13°  to  19°  inwards  towards  Neptune  as  measured  from  the  tangent 
to  the  original  circular  orbit.  Pluto,  given  these  conditions  Immediately 
starts  to  exhibit  the  characteristics  of  resonant  oscillation  and  stability 
mentioned  before. 

First,  since  the  escape  must  take  place  in  Neptune’s  orbit  the  escape 
point  is  both  the  node,  the  point  of  conjunction^ and  near  the  perihelion 
of  Pluto's  new  orbit.  It  was  shown  that  if  the  conjunction  occurs  any- 
where other  than  the  the  aphelion  the  resonance  will  oegin  to  oscillate 
about  the  aphelion.  Fbr  the  case  shown  in  Figure  13  it  can  be  seen  that 
the  resonance  is  initially  a 3/2  resonance.  The  oscillations  are  very 
apparent.  Because  they  start  near  the  perihelion  the  amplitude  is  large. 
Figures  13,  14,  and  15  show  the  period  from  escape  at  t ■ 0 out  to  135 
orbits  of  Neptune.  It  can  been  seen  that  the  oscillations  are  slightly 
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FIGURE  13.  Oscillating  3/2  Resonance  from  0 to  45  orbits 

after  Pluto’s  escape  from  Neptune 


17  DEG  DBASIC  RUN  H/U,S,J  J-123 

FIGURE  14,  Pluto' 8 3/2  resonance  pattern  froa  45  to  90  orbits  of 
Neptune  ( 15,000  yrs.)  after  the  escape. 
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damped  and  pluto  at  perihelion  doeB  not  quite  return  to  the  node.  Thus, 
at  least  for  a short  period  the  resonance  Is  stable. 

One  of  the  primary  characteristics  of  an  orbital  resonance  discussed 
Is  the  fact  that  once  a resonance  Is  established  all  traces  of  any  Initial 
conditions  are  lost.  Integrating  a resonance  backwards  and  forwards  In 
time  will  not  produce  any  change  In  the  overall  orbital  configuration, 

A test  of  the  escape  formed  Neptune-Pluto  3/2  resonance  would  be  to  inte- 
grate from  a point  hack  towards  the  moment  of  escape.  If  the  two  planets 
continue  to  orbit  the  Sun  in  resonance  past  the  time  of  escape  then  It 
will  strongly  suggest  that  a resonance  has  successfully  been  formed  and 
will  show  that  no  information  as  to  a resonances  Initial  conditions  can 
be  derived  from  that  resonance.  This  process  Is  shown  in  Figures  16, 

17,  1^,  and  19  which  show  Integrations  from  a point  120  orbits  after 
the  escape  back  to  a point  -60  orbits  before  the  escape.  Comparing 
these  with  Figures  13  through  15  we  see  that  the  resonance  pattern  evolves 
when  integrating  forward  in  time  and  then  when  we  reverse  the  process 
we  can  never  return  to  the  escape.  This  strongly  suggests  that  at  least 
lnltla'ly  Pluto  can  escape  from  Neptune,  onter  a 3/2  resonance  which 
oscillates  about  the  aphelion  and  contains  no  Information  concerning 
Its  Initial  conditions. 

The  orbital  characteristics  of  Pluto  after  escape  are  Interesting 
also.  The  orbital  inclination  was  taken  to  be  defined  as  17°.  The 
eccentricity  varies  sllgHUy  with  each  orbit  under  the  perturbations 
but  Just  from  the  escape  conditions  necessary  to  produce  the  resonance 
the  orbit  has  an  eccentricity  of  e ■ 0.24.  The  compares  with  the  measured 
value  of  e - 0.249.  The  period  and  the  semi -major  axis  are  governed  by 
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Kepler’s  Third  Ian  and  with  the  conditions  necessary  for  resonance  we 
arrive  at  valuea  for  period  and  semi -major  axis  that  are  close  to  the 
present  measured  values.  See  ftgure  20. 

Thus  it  is  possible  for  Pluto  to  escape  from  Neptune  and  go  into 
a resonant,  orbit  with  Neptune  considering  only  the  Sun,  Neptune  and 
Pluto,  Gi *en  the  proper  initial  conditions  that  resonance  will  be  a 
3/2  type  and  will  have  nost  of  the  characteristics  that  are  measured 
today. 

The  next  step  is  to  add  in  the  other  planets  in  the  outer  solar 
aysten  in  order  to  determine  what  effect  they  have  on  the  stability  of 
the  resonance  and  its  evolution  into  Its  modern  configuration. 
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Effects  of  Uranus  on  the  3/2  Resonance 

B y Integrating  the  notion  of  Pluto  under  the  influence  of  the  Sun, 
with  Neptune  as  the  only  perturbing  force,  we  have  seen  that  a 3/2  res- 
onance is  at  least  possible.  Due,  however,  to  the  large  amplitude  of 
the  oscillations  and  the  proximity  of  the  perihelion  and  the  point  of 
conjunction  to  the  orbital  node  the  resonance  lasts  only  a few  complete 
oscillations  before  it  Is  perturbed  into  a slightly  different  non-res- 
onant orbit.  Figure  21  shows  this  sort  of  perturbation.  Something 
is  needed  to  stablize  the  resonant  oscillations. 

The  Sun,  Neptune  and  Pluto  are  a rather  crude  approximation  of 
the  outer  solar  system.  Upon  close  examination  Uranus,  the  next  planet 
inward  from  Neptune,  is  very  close  to  being  in  a 2/1  resonance  with 
Neptune  and  consequently  would  be  in  a 3/1  resonance  with  Pluto.  Adding 
Uranus  to  the  model  solar  system  has  an  immediately  visible  effect  on 
the  evolution  of  the  resonance.  Whereas  before,  Pluto’s  conjunction 
with  Neptune  oscillated  nearly  180°  with  very  little  damping,  now  with 
Uranus  present  the  amplitude  of  the  oscillations  slowly  decreases.  This 
Cart  be  seen  by  comparing  Figure  22  with  Figure  15.  The  "turn-around" 
point  of  the  oscillations  is  clearly  displaced  away  from  the  node  on 
the  x-axis.  The  stable  point  is  on  the  y-axis  representing  the  aphelion 
conjunction.  Thus  the  yf\  resonance  between  Pluto  and  Uranus  has  the 
combined  effect  of  damping  the  oscillations  so  that  the  resonance  begins 
to  take  on  some  of  its  present  day  characteristics  and  in  doing  so  Uranus 
prevents  Pluto  from  ever  coming  close  enough  to  Neptune  to  be  perturbed 


FIGURE  21.  Diagram  showing  a 3/2  resonance  being  perturbed 

out  of  the  resonance  and  into  a non-resonant  orbit 
Rotating  Coordinate  Frame 
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| out  of  the  orbital  resonance.  This  Is  the  situation  today  when  Pluto 

never  comes  closer  than  16  A.U.  to  Neptune.^  We  now  have  an  escape 
! formed  3/2  resonance  which  is  independent  of  the  initial  start  that  is 

g stable  with  Uranus  and  Neptune  for  approxiamately  1340  orbits  of  Neptune 

or  the  equivalent  of  225»000  years.  After  that  the  resonance  begins 
| to  show  signs  of  decay  and  the  amplitude  of  the  oscillations  Increases 

until  Pluto  is  perturbed  out  of  the  resonance.  Long  term  perturbations 
J app®ar  to  be  needed  to  give  the  orbit  additional  stability, 
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toiw  Term  Effseta  of  Saturn  and  Juplitr 

Even  a vary  alapla  look  at  the  solar  systsa  sufflesa  to  revsal  that 
ths  major  henstituents  ara  the  Sun  and  Jupltsr,  Thus,  In  any  study  In- 
volving perturbations  among  ths  planets  we  must  include  the  influence 
of  Jupiter  and  Saturn,  Although  Jupiter  and  Saturn  are  a much  greater 
distance  from  the  Pluto-Neptune  system  they  more  than  make  up  for  this 
with  their  large  masses.  Because  of  their  shorter  periods  of  revolution 
their  influence  is  mainly  over  a longer  period.  This  is  exactly  what 
is  needed  to  add  to  the  stability  of  the  resonance. 

The  first  noticeable  effect  of  the  presence  of  Jupiter  and  Saturn 
is  that  Jupiter  with  its  extremely  large  mass  can,  if  positioned  correctly, 
perturb  Pluto  out  of  the  reaonanoe.  This  occurs  only  during  the  initial 
stages  of  the  resonanoe  formation.  If  Neptune  and  Jupiter  are  in  conjunc- 
tion at  a time  when  Pluto  is  also  near  conjunction  and  perihelion  passage 
ths  combined  perturbing  forces  are  enough  to  perturb  Pluto  out  of  the 
resonanoe.  Jupiter  must  be  positioned  so  that  it  will  not  be  in  such 
a position  until  the  resonance  has  had  a chance  to  shift  enough  to 
prevent  Pluto  from  approaching  closely  to  Neptune.  If  Jupiter  is  placed 
so  that  after  Pluto's  escape  from  the  Influence  of  Neptune,  Jupiter  and 
Pluto  are  in  conjunction  when  Pluto  is  at  or  reasonably  near  its  greatest 
distanoe  relative  to  Neptune,  then  Jupiter  will  not  perturb  the  resonance 
and  the  long  term  perturbations  due  to  Saturn  and  Jupiter  can  be  fully 
evaluated. 
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With  the  Addition  of  Jupiter  and  Saturn  we  now  have  a complete 
computer  model  of  the  outer  solkr  system  and  a very  good  approximation 
of  the  entire  solar  system  due  to  the  small  relative  masses  of  the  inner 
terrestrial  planets.  We  can  now  study  the  evolution  of  the  Pluto - 
Neptune  resonance  in  detail. 

It  was  found  that  if  the  escape  of  Pluto  occurs  in  such  a manner 
that  when  Pluto  is  initially  at  its  fartheet  relative  distance  from 
Neptune  it  is  near  conjunction  with  Jupiter  then  a stable  3/2  resonance 
can  develop.  Pluto’s  new  orbit  has  an  eccentricity  roughly  equal  to 
0.24.  The  resonance  oscillates  with  a very  large  amplitude  about  the 
aphelion  conjunction  point.  This  pattern  was  found  to  be  constant  for 
nearly  7000  orbits  of  Neptune  which  is  the  equivalent  of  1,176,000  years. 
There  are  no  signs  that  the  resonance  can  be  perturbed  out  of  this  con- 
figuration, The  oscillations  have  been  damped  enough  so  that  they  no 
longer  come  near  Neptune.  In  addition,  as  was  mentioned  before,  reversing 
the  integration  produces  no  changes  in  the  resonance.  Thus,  with  all 
the  planets  present  we  are  able  to  form  a true  resonance  between  Pluto 
and  Neptune. 

One  additional  feature  of  the  resonance  was  identified  as  a result 
of  long  term  integrations.  examining  the  velocity  in  the  z-direction 
at  the  moment  of  perihelion  passage  enables  us  to  locate  the  perihelion. 

At  escape  the  perihelion  z-velocity  was  near  its  maximum  value  indicating 
that  the  perihelion  lay  near  the  node  between  Pluto's  and  Neptune's 
orbits.  As  the  integrations  progressed  the  z-velocity  at  perihelion 
decreased  to  zero  and  then  began  to  increase  in  the  negative  direction. 
This  indicated  that  the  perihelion  was  precessing  away  from  the  node 
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and  up  away  from  the  plane  of  the  ecliptic.  Current  studies  have  shown 

that  the  perihelion  oscillates  with  a period  of  four  Million  years  about 
o 7 

a point  90  away  from  the  orbital  node.  Upon  carrying  out  the  integration 
to  slightly  over  one  Billion  years  we  found  that  the  perihelion  was 
circulating  Instead  of  oscillating.  There  have  been  studies  that  show 
am  circulating  perihelion  but  they  do  not  include  Jupiter  and  Saturn 

Q 

in  theii  analysis  of  the  Pluto-Neptune  system,  i'he  fact  that  the  peri- 
helion circulates  rather  than  oscillates  does  not  invalidate  the  other 
results.  It  is  certainly  possible  that  continued  integrations  or  more 
refined  calculations  of  the  perturbations  due  to  Jupiter  and  Saturn  might 
change  the  circulating  perihelion  to  its  present  oscillating  form. 

Ons  -Million  years  is  a long  time  on  man’s  sclae  but  a mere  flicker 
astronomically. 


DBASIC  RUN  U/U.S.J  J-123  1010-1033 

FIGURE  23.  Pluto-Ifeptune  3/2  resonance  1055  orbits  after  escape 


5495-5540  17  DEG  DBmSIC  RUN  W/U.S,J  J-123  0=3350 

FIGURE  25.  Pluto  and  Neptune  3/2  resonance  924,000  years 

after  Pluto  escape  fro*  Neptune 
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has  .been  shown  that  it  is  possible  that  if  Pluto  and  Triton  had 
a close  encounter  .resulting  in  Pluto  attaining  escape  velocity  directed 
approximately  17°  inward  from  its  previous  orbit  then  a stable  3/2 
oscillating  resq^noe  will  evolve.  Pluto  will  orbit  the  Sun  in  an  orbit 
that  has  all  the  characteristics  of  its  present  day  orbit.  The  resonance 
has  been  found  to  be  stable  for  over  one  millior.  years  with  no  evidence 
to  suggest  that  it  will  not  continue  to  remain  stable.  Pluto's  perihelion 
was  found  to  circulate  rather  than  oscillate  but  this  is  explainable  as 
a result  of  the  perturbations  from  Jupiter  and  Saturn. 

The  objection  may  arise  that  we  have  placed  too  many  restrictions 
on  the  solution  and  that  the  solution  may  have  been  forced.  This  ignores 
the  fact  that  any  resonance  in  itself  is  a special  case  and  a 3/2  resonance 
is  a special  ease  of  that.  Many  initial  conditions  were  found  that  allowed 
Pluto  to  escape  and  orbit  the  Sun  without  entering  a 3/2  resonance.  But 
Pluto  is  in  a 3/2  resonance  today  so  we  had  to  look  for  the  particular 
conditions  that  gave  us  a 3/2  resonance. 

This  project  was  carried  out  for  the  purpose  of  determining  whether 
it  was  possible  for  Pluto  to  have  escaped  from  Neptune  and  evolve  into 
3/2  resonance.  Such  a possiblity  was  found  and  in  doing  so  we  have 
eliminated  the  major  objection  to  the  moon  hypothesis  as  to  Pluto's  origin. 
That  objection  was  that  since  Neptune  and  Pluto  are  in  resonance  they 
would  never  come  close  enough  to  one  another  for  Pluto  to  have  been  a 
moon.  We  have  shown  that  an  escape  will  produce  a resonance  and  once  the 


resonance  is  established  all  the  information  regarding  the  initial  start 


is  lost.  The  evolution  of  the  resonance  keeps  Pluto  and  Neptune  separated 


and  has  all  the  modern  characteristics  that  have  been  observed  in  the 


system 
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Appendix  I 

A Method  for  Determining  Resonant  Stabilities 
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One  of  the  interesting  side  results  of  this  investigation  has  been 
the  discovery  that  certain  resonances  are  more  stable  against  perturbations 
in  the  oscillations  than  others,  this  also  has  a definite  bearing  on  the 
total  number  of  resonant  combinations  that  can  exist.  We  defined  a 
resonance  as  any  system  where  the  periods  of  the  bodies  were  the  ratio 
of  two  small  whole  numbers.  Taking  this  to  be  any  number  less  than  ten 
and  eliminating  multiples  and  so-called  mirror  resonancesi  that  is,  a 
3/1  resonance  is  the  same  as  a 1/3  resonance  except  viewed  from  the  per- 
turbing planet  instead  of  from  the  perturbed  planet  as  usual,  we  are 
still  left  with  2?  different  combinations. 

While  each  resonance  is  stable  in  its  non-oscillating  formt  if  it 
is  perturbed  or  is  oscillating  a majority  of  these  patterns  get  perturbed 
into  another  orbit,  A possible  explanation  for  this  can  be  found  by 
going  beck  to  the  basic  mechanism  that  maintains  a resonance,  the  oscill- 
ations of  the  conjunction  point.  These  oscillations  show  up  as  small 
variations  in  the  relative  period  of  the  outer  planet  relative  to  the 
inner.  Fbr  a 3/2  resonanoe  the  period  between  perihelion  passages  is 
1.50  orbits.  Similarly,  for  a 4/3  the  period  is  1.33  and  for  an  8/7 
the  pe-iod  is  1,14,  or  in  general  the  period  is  just  the  numerical  value 
of  the  resonance  ratio.  The  oscillations  slow  or  speed  up  these  periodst 
typically  for  a 3/2  resonance  the  values  range  from  1.47  to  1.53.  Right 
away  we  see  that  a resonance  must  have  a certain  amount  of  "room”  to 
oscillate  in  for  it  to  be  distinguishable  as  an  oscillating  resonance. 

If  there  was  another  resonance  within  the  range  of  the  oscillations  there 
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would  be  some  sort  of  interference  between  the  two. 

This  is  the  situation  for  an  8/7  resonance.  It  is  stable  in  the 
non-oscillating  case  but  if  it  oscillates  it  interference  with  the  9/8 
at  1.13  on  the  low  side  and  the  7/6  at  1,17  on  the  high  side.  In  addition, 
the  8/7  resonance  has  a loop  that  falls  near  the  perturbing  planet  so 
that  any  movement  brings  this  loop  and  thus  the  planet  closer  to  the 
perturber  and  consequently  increases  the  disruptive  force  of  the  planet. 

The  amount  of  relative  "room"  between  resonances  can  be  seen  if  we  plot 
the  relative  periods  on  a number  line  and  then  look  for  gaps  in  which 
a resonance  could  oscillate.  As  it  turns  out  the  gaps  occur  around  the 
resonances  3/2,  5/3 * 4/3*  2/1,  3/1,  5/2,  and  slightly  less  stable  ones 
around  5/4,  7/3,  and  8/3  resonances.  These  are  the  same  resonances  that 
we  see  throughout  the  solar  system.  Virtually  none  of  the  other  resonances 
are  found. 

Looking  for  the  relative  amount  of  "room"  around  a resonance  is 
only  a convenient  way  of  finding  resonances.  There  is  more  likely 
a basic  explanation  of  the  resonances  but  is  interesting  that  this 


method  works  at  all 
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The  following  computer  program  is  a model  of  the  Neptune  satellite 


system  while  it  is  orbiting  the  Sun.  It  was  used  to  study  the  initial 
escape  of  Pluto  from  orbit  around  Neptune.  The  program  is  made  up  of 
four  partsj  a force  function  definition  section,  an  initialization  of 
variables  section,  an  integration  section  using  a Runge-Kutta  numerical 
integration  scheme,  and  finally,  an  output  section. 


'1 

. 


j ; 

ii 


meptune 
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2UREM  T IE  FOLLOWING  STATEMENTS  SF.T  UP  OUTPUT  FILES 
21FILE  ‘f  1 : " YVX"  'Y  \/S  X 
22FILE  #2:  "XVT"  'X  VS  T 
23  FILE'»3:"3!>" 

25SCR ATCU  'M 
26SCRATCU  #2 
27  SCRATCHES 
99’ 

10URE-1  THE  FOLLOW  I MG  STATEMENTS  ARE  ACCELERATIONS 
105LET  Pis  12. 566370 61 4 4 
110  LET  P2=39.47S4176e44 
120DEF  F N A ( A , B , C , D , E ) 

1 2 1 LET  FMAs  -764966 . 42* ( A+1 1 O83 -765 ) /( ( ( A+ 1 19*3 . 755) ~2+B~2) * ( 1 . 6) ) 

122  LET  FMA  = FMA-P2*A*(A*A  + R*B)“(-1.6)+4.46337ME-7*(A+11933.766) 

123  LET  F‘1A  = FNA+.U01334673*C 

124 

129FMENI) 

1 3 e D E F F N B ( A , B , C , D , E ) 

131  LET  FMB--764966 . 4?*B*  ( ( A + 1 1 9S3 . 75S ) “2+B*  B)  ~ ( - 1 .S)-P2«B*(  A*A+H»B)  “( 

132  LET  FN3=FNB+4 . 463377  3F-7*  3- .00 1 334673*C 
1 39FMEMD 

14c  :)EF  FMC  ( A , 3 , C , 0 , E) 

141  LET  FMCs-P2»E*(  A"2+B"2+E‘*2)  **  ( - 1 .6) 

142  IF  MrO  THEM  149 

143  LET  Q1 s A- 1-( 1-D) * ( COS( 6 .2331 35307*T0) -1 ) 

144  LET  Q2- B-( 1 -D) * SIN ( 6 . 283 1 86307*T0 ) 

146  LEI  FNC--FNC-P2*  1*E*  ( 0 1 ''2+02~?+E~  2)  ~ ( - 1 .6) 

149  r NE  JO 
199’ 

2UUREM  THE  FOLLOWING  ARE  INITIALIZING  STEPS 

2e1PRlNT"IMPUT  RATIO  OF  PERTURBING  TO  PRIMARY  MASS.  TF  THERE  IS  NO" 
2e2PRI NT" PERTURBING  MASS,  INPUT  0" 

20 31 M PUT  M 

2o6PRliNT"  INPUT  0,1  FOR  NON-ROTATING  OR  ROTATING  COORDINATES" 

207  INPUT  1) 

209PR1 NT" INPUT  INITIAL  X,  VX" 

2101 NPUT  X,U 

2 1 5PKI NT" 1 NPUT  INITIAL  Y,VY" 

220  IN  PUT  Y,V 

222  PRINT  "INPUT  INITIAL  7.,V 

223  INPUT  Z,W 

226PI1 1 NT"  INPUT  STARTING  TIME,  END  TIME,  ft  STEPS  .OUTPUTSTEP  INTERVAL  " 

2 3c  INPUT  To , T 1 , N , M 1 

233LET  H-Cri-TO/M 

246  LET  R SQ W ( X * 2+ Y ~ 2+Z~  2 ) 

246  LEI  T2a ( 1 So* ATM ( Y/X ) ) /3 . 1 4 169 

247  IF  X>0  THEM  260 
24d  LET  T2-T2+180 

260  PRINT  ^ 3 : X , " , " , Y , " , " , Z 
266PR1NT  112:  To  , X 
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NEPTUMF.  (continued) 


IIS  PAGES  IS  BUST  QUALITY  PRACJICAiL* 
ISOM  COPY  FURNISHES)  10  DD.O 


399' 

400 REM  DIE 
401  FOR  Is 
406 R KM  X 1 
4 K ...  ET  K1  = 
4 1 5LET  !<2r 
417  LET  '(5 
4 2oL F.T  K3- 
4251 ET  <4- 
427  LET  K6 
4 30 LET  Ills 
4 V > i.  E T VI- 
4 37  L ET  •■'1 1 
4401  ET  XU 
443LET  YU- 
447  LET  Z1 
4 50 LET  02- 
4 56LF.T  V2-- 
457  LET  W2 
4 6 0 L E T X2- 
465LET  Y 2- 
467  LET  12 

4 7 0 L L T K 1 - 
475LET  K2r 
477  LET  K5 
4'k'LET  !<3- 
4 36 LET  K4- 
437  LET  K5 
4Q0LET  UU 
495L.ET  Vis 
497  LET  VI 
5 00 LET  XU 
5051  ET  YU 
30 7 LET  7.1 

5 10l.  c.T  U2- 
316LET  V2s 
617  LET  W 2 
520LET  X2s 
525LET  Y2s 
527  LET  72 
5 30 LET  K 1 -• 
535LETK2 - ! ! 
537  LEI  K5 
540LET  K 3- 
5 4 5 L E T 1(4  r 
547  L E T L.  6 
550LET  tJH- 
555LET  VU 
557  LET  >;  1 
5 6 0 L E T XI- 


f OL LOW  1 90  STFU3  PERFORM  PIE  1 9T  EG  RATIO'! 

1 TO  II 

GOES  a'ITH  '.),  K2  WITM  V,  K3  WIT'!  Y,  K4  WTTM  Y,'<6  WITH  V, 

M*F)IA(X,Y,V,D,7) 

il* FN 9(  X , Y , U , D , Z ) 

s ! 1 * F ’ I C ( X , Y , V! , D , 7. ) 

!)*U 
il»V 
-M»  V 
l)+<1/6 

V + K 2/6 

- i+<S/6 

x + -;  3/4 

Y+  (4/(1 
sZ  + '(6/  6 
U+X1/2 
V+K2/2 

s W + '<  6 / 2 

X+K3/2 
Y + '(  4 / 7 
sZ  + K5/2 

i|*FNA(X7,Y7,V2,D,7.2) 

H * F 9 B ( X 7 , Y 2 , U 2 , D , 7. 2 ) 
sM*F:iC ( X2  , Y2 , V2,D,7.4>) 

H»U2 
H*V2 
sG*  72 
l)  1+'<  1 / 3 
vi  +'<7/ 3 
s V1+K6/3 
X 1 +X  3/ 3 
Y1+X4/3 
-Z1+K5/3 
U+X 1/2 
V+K2/2 
s 7 + K5/2 
X + '<  3 / 7 

Y + ’<  4 / 2 
-Z  + (6/2 

H*F'IA(X?,Y2,V2,D,7.7) 

*FNn(X2,Y2,U7,n,Z2) 

- T * F M C ( X 7 , Y 7 , 7 7 , )) , 7 '■ ) 
ll*U2 

:i*V2 
- 1 1 * >/  2 
U1+K1/3 
VI +'<2/ 3 

- VUl(5/3 
XU'O/1 


X6VITF"'7 
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IBIS  ?AQS  IS  B«ST  TWmCaa* 
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NEPTUNE  (continued) 

565LET  YUYWX4/3 

567  LET  Z 1 - Z 1 +K6/3 

57eLET  U2-U+K1 

575LF.T  V 2 - V + K 2 

577  LET  W2-W+K5 

5SCLET  X2sX+K3 

535LET  Y 2 - Y + K 4 

537  LET  Z2&Z4-K6 

59c LET  K1-  1»FNA(X2,Y2,V?,D,Z?) 

595LET  K2-ll*F‘JB(X2,Y2,l)2,D,Z2) 

597  LET  K5sil*FNC(  X2 , Y2 , V2 , D , Z2) 

6c c'LET  K3sU»l)2 
5 c 5 L E T K4s:H#  V2 
6 >.7  LET  K6-M»>f2 
61CLET  U- U 1 +K 1 /6 
615LET  VsV1+X?/6 
517  LET  ./s V1+K5/6 
62CLET  XvX1+K3/6 
625LET  Y-YUK4/5 
627  LET  Z-ZUK6/6 
63c LET  TcsTO+M 

6351F  INT(I/M1)-(I/M1)Oc  THEM  665 
640  LET  R-SQR(X<'2+Y‘'2+Z<‘?) 

6*4 1 LET  T2s(  1’U'*ATN(Y/X)  )/3.  14159 

642  IF  XX'  THEM  645 

64  3 LET  T2-T2+1.3C 

645PHIMT  n:X,Y 

65cPRi:iT  '/ 2 : T c , X 

665  PKI’II  >* 3 : X ; " , " ; Y ; " , " ; 7. 

67C  ME  XT  I 

1 14c?KnT"FIIIAL  T I ‘ 1 F.  - " ; Tc 

1 145PR1  TP'FIMAL  X-  ";X;"VXr  " ; IJ 

1147  PHI  NT" FINAL  Y. " ; Y ; " VY- " ; V 

115C  PRINT  "FINAL  7.  = ";  7.;"  VZ-";V 

HCOOHESET  *1 

He  1 0 RESET  It 2 

HC3C  RESET  II 3 

90UCPRINT 

9999EN0 


The  following  program  Is  the  complete  model  of  the  outer  solar 
system  tontainlng  all  the  plants  from  Jupiter  outwards.  The  program 
can  be  run  In  both  a Cartesian  fixed  coordinate  frame  or  In  a rotating 
relative  coordinate  frame.  The  name  LONCZUSJ  is  derived  from  the  fact 
that  the  program  was  designed  to  study  the  3/2  Pluto  - Neptune  resonance 
over  long  time  scales  and  while  doing  so  would  print  out  the  ^-position 
and  ^-velocity  of  Pluto  each  time  Pluto  was  at  perihelion  so  that  we 
could  study  the  precession  of  the  nerihelion.  In  addition,  the  program 
contains  the  planets  Uranus,  Saturn,  and  Jupiter.  Each  planet  Is  defined 
by  its  position  in  space  In  each  of  the  three  coordinates  using  only 
three  lines  in  the  program.  FY>r  example,  Sa tarn's  position  Is  computed 
by  lines  293-2°^,  393-39^,  end  493-404.  Saturn's  perturbation  on  Pluto 
Is  then  computed  and  processed  at  each  step  of  the  Runge-Kutta  Integration 
sequence  with  Is  almost  Identical  to  that  in  the  program  NEPTUNE. 
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TiOWGCTBJ 

1 DIM  P(IOO) 

2 SETDIG1TS  7 

20  LET  P7  = 35* . 1 4SS70U7 

100REM  THE  FOLLOWING  STATEMENTS  SFT  UP  OUTPUT  FILES 
110  FILE  1 1 : "INITI AL2" 

120  FILE  #2: "DATA V 
130  FILE  03:"XVX1" 

1 35FILE'<4 : "XY3" 

140’ 

150RE3  THE  FOLLOWING  STATEMENTS  ARE  ACCELERATIONS 

160  LET  P 3s. 6 37 86 

170  LET  P4=12. 32461 

ISO  LET  P5-4.360SF.-5 

185  LET  PO*. 1729238 

19ULET  PI-12. 5563706144 

195  LET  P6-.3148958R7 

200  LET  P2- 39. 4784  176044 

202  LET  P( 10)sR7. 28965 

205  LET  PS-2.S5SE-4 

207  LET  P(  1 1 )=').547966E-4 

210DEF  FNA(A,D,C,D,E) 

220LET  FNAsP2*(D-(  A''2+D'*2  + E'*2)  " ( - 1 . 5 ) ) * A+P 1 »C*D 
2 30 IF  HeO  THEN  300 

240LET  Q U- A- 1-( 1-D)»( COS( 6.2831 85307»TC>-1) 

250LET  02 sB-(1-D)*Sl'l (6. 28318530 7*T0) 

260LET  FMAsFNA-P2»Q1*M'»(Qr2+Q2“2+E“2)“(-1 .5) 

270  LET  Q3*A-P3»(COS((P4-D*6.283185307)»T0)) 

2 8o  LET  04 i- 8-P3*SI t!((PN-D* 6. 28318530 7)*T0) 

29o  LET  F“J AsFU A-P?*03* P6* ( 03*2+04*2+E'*2 ) * ( - 1 . 6 ) 

293  LET  05* A-P6* ( COS( ( P7-D*6 . 2831 R5307 ) •TO) ) 

294  LET  0S-!UP6»SIN((P7-D*6.2331853O7)*To) 

295  LET  FNA a FN A-P2*05* PS* ( f)6*2^08''2+E“2 ) * ( - 1 . 5 ) 

296  LET  ')7*A-P9»(COS((P(  1 o ) -D»6 . 25  3 1 85307 ) •TO- 1 23 . 2285 ) ) 

297  LET  ()8s1ur»9»Sl'l(  (P(  1 0 ) -0*6 . 28 3 1 85307 > *T0  -123.2285) 

298  L ET  FNA«FNA-P?*  97*P(  1 1 ) » ( Q7"2+0S‘*?+E‘*2>  * ( - 1 . 5 ) 
300'FME.Nf) 

310DEF  FNH(  A , H,C , I) , F. ) 

320L.ET  FNI1*:P2»(D-(A*2+H*2+F>‘*2)*(-1.5))»R-P1»C»D 
3301F  MsO  THEN  400 

340LET  QUA-1-(  1-D)»(COS(6.283185307»TO)-1) 

3 50 LET  Q2sB-(  1-D)»SIN(6.2R3185307*T0) 

360L.ET  FNHsFNB-P2#Q2*M* ( 01  ‘*24.r)?‘* 2+F.‘'2)*(-1 .5) 

370  LET  03  = A-P3#(COS((pt»-D»6.283185307)»T0)) 

380  LET  04i.n-P3»SIN(  ( P4-D»6.?.33185307>#T0) 

390  LET  FMru-FNH-P2»94»P5»(93‘‘2*04"24.E~?)'‘(-1 .5) 

393  LET  05  s A-P6  • ( CO.S(  ( P7-D*S  , 28  3 1 85  307 ) *T0 ) ) 

394  LET  ()6sB-PS»SIN((P7-n»6.?83lS5307)»T0) 

395LET  FMHi.FNU-P2*'')6# PR » ( 05*2+OS‘'?+E'* ? ) “ ( - 1 . 5 ) 

396  LET  07sA-P9*(COS( (P( 1o)-D*6. 283 185307 )*T0- 123. 2285)) 
ovfl 
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398  LET  FN9rFNB-P2*Q8*P(  1 1 ) * ( 07"?+Q3~2+E~2) * ( _ 1 .5)  ~ 

4 OOF MEND 

41o  DEF  FNC(A,n,C,D,E) 

42 0 LET  FNC=-P2*F.»(  A“2+n*2+E“2)“(-1 .5)  ~ 

430  IF  MsO  THEM  500 

440  LET  0 1 - A - 1 - ( 1 - D ) * ( COS  ( 6 . 8 3 1 8 5 3 07  * TO ) - 1 ) 

450  LET  Q2  = !U(1-D)*SlN'(6.?331353e7vro) 

4 So  LET  FNC  = FriC-P?*M»E»(or'?«-02~?+F.~?)~(-1 .5) 

4 ?0  LET  Q3rA-P3»(COS((P4-D»4.2a?135307)*T0)> 

4 80  LET  Q4=B-P?»SIN(  (P4-n4*5.2831S5307)*TO) 

49o  LET  FNCsFMC-P?*E*P5*(03‘*?+'54"24.F.‘*2)“(-1 .5) 

493  LET  Q5  = A-P64 ( COS ( ( P7-9*6 . 283 1 85  307 ) *T0 ) ) 

404  LET  Q6s.-P.-P6*SIM(  (P7-9*6.?.331S5307)*T0) 

4)5  L E T F N C - F '1 C - P 2 * S * P 8 4 ( 0 1 5 ** 2 + 9 1 6 “ ? ♦ F **  2 ) **  ( - 1 . 5 ) 

4 9 5 LET  07= A-P9»(C03( ( P( 1 0 ) -0*5 . 28  3 1 85  307 ) *T0- 1 23 . 2235 ) ) 

497  LET  Q3  = !3-P9*SlN((P(1O)-D’»6.233185307)*TO  -123.2235) 

493  LET  FHC  = FMC-P2*  F.*  ? ( 1 1 ) * ( 97  ''2+03‘*2+F.'*2)''(-1 .5) 

5 00  FNEMD 

510'  ~ 

520REM  T!1E  FOLLOWING  ARE  INITIALIZING  ST^.PS 

5 3o  PR  I NT”  INPUT  RATIO  OF  PERTURBING  TO  PRIMARY  MASS.  IF  TMF.RF  IS  NO" 

540 PRINT" PERTURBING  MASS,  INPUT  0" 

5 50  IN  PUT  M 

560 PR  I NT" INPUT  0,1  FOR  NON-ROTATING  OR  ROTATING  COOPDINATFS" 

57 0 INPUT  D 

580 PR  1 NT" I N PUT  INITIAL  X,  VX" 

5 90  IN  PUT  X,IJ 

SooPR  I NT" 1 M PUT  INITIAL  Y,VY" 

610INPUT  Y , V 

420  PRINT  "INPUT  INITIAL  7,V7." 

630  INPUT  Z ,W 

640 P R 1 NT"  1 N PUT  STARTING  TIME,  END  TIME,  « STEPS , OUT  PHTSTF.P  INTERVAL " 

6501 N PUT  TO ,T1 , H , N1 
660LET  !U(ri-TO)/N 
670  PRINT'*?:  X,U 
6 So  PRINT#? : Y , V 
69o  PRINT#2: Z , W 
7oo  • 

710REM  THE  FOLLOWING  STEPS  PERFORM  THE  INTEGRATION  ^ 

720  FOR  1=1  TO  M 

730REM  K1  GOES  7ITM  U,  '<2  WITH  V,  K3  WITH  v,  K4  WITH  Y,K5  WITH  ’■'.WTTM  7 
740L.ET  KUH»FNA(X,Y,V,D,Z) 

750LET  :C2  = H*FNn(X,  Y,U,D,7) 

760  LET  K5=H*FMC(X,Y,W,n,Z) 

770LET  K 3 - ' I * U 
780LET  K4=H»V 
790  LET  K6=U*W 

300LET  U 1 =U+< 1 /6  ^ 

310LET  VUV+K2/6 
320  LET  WU-W+K5/6 
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LOKZUBJ  (continue!) 

330LET  X 1 ^X+K  3^6 

840LET  YUY+K4/6 

850  LET  7 1 sZ+K6/6 

36oLET  U2iU  + '<  1 /2 

870LET  V2i-V+X2/2 

380  LET  W2sW+'<5/2 

390LET  X2-X+X3/2 

900LET  Y2- Y+K4/2 

910  LET  Z2sZ+K5/2 

920LET  ’<U-'1*FNA(X2,Y2,V2,n,Z?) 

9 3 0 L E T K2=!1*FN3(X2,Y2,U2,D,Z2) 

940  LET  :<5-H*F:lC(X2,Y2,W2,Df  7.2) 

950LET  f< 3 - H » U 2 

96CLET  K4sH*V2 

970  LET  !<6s:H»'72  * 

980LET  UUUUK1/3 

9901. ET  V1-VHK2/3 

1000  LET  WUWUK5/3 

1010LET  X1-X1+K3/3 

1020LET  Y1-Y1+K4/3 

U30  LET  Z1-Z1+K6/3 

104CLET  U2-U  + X 1/2 

1O50LLT  V2-V+K2/2 

1060  LET  «.'2=W+K5/2 

1070LET  X2-X+K3/2 

I080LET  Y2-Y+;<4/2 

1090  LET  Z2-7.+K6/2 

1 100LET  K1-H»FMA(X2, Y2,V2,n,7.2) 

1 1 1 0 L E T K 2 - H * F M B ( X 2 , Y 2 , U 2 , D , 7. 2 ) 
1120  LET  K5-M»FNC(X2,Y2,W2,n,72) 
1130LET  K3-n»lJ2* 

1140LET  K4 -H»V2 

1150  LET  L6sllM2 

1 1 60 LET  UUU1+’<1/3 

1170LET  V 1 -V 1 +K2/3 

1130  LET  VUW1+K5/3 

1190LET  XU-X1+K3/3 

1200LET  Y 1- Y 1 +!<4/3 

1210  LET  7 UZ  1+1(6/ 3 

1220LET  U2-U+K1 

1230LET  V 2 ::  V ♦ K 2 

1240  LET  /^2iW+K5 

1250LET  X 2 i- X 4- K 3 

1260 LET  Y2i: Y + K4 

1270  LET  Z2iZ+K6 

1 230L ET  '< 1 sill »FN  A ( X2 , Y2 , V2 , 1) , 72 ) 

1290LET  K2s»l»F?in(X2,Y2,lJ2,D,Z2) 

1 300  LET  K5ril»FNC(X2,Y2,U2, 0,7.2) 
1 3H LET  <3-H»U2 
1320LET  K4-I1  *V 2 
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LONGZUSJ  (continued) 


1330  LET  K6s»l»W? 

1 340LET  U;Ul+l(1/() 

1330LET  V-V1  + K2/6 

1360  LET  W-N1+K5/6 

1 370 LET  XsXUK3/6 

1 3'K'LET  Y-Y1+K4/6 

1 390  LET  Z-Z1+’<6/6 

1400LET  TO-TO+H 

14  10  LET  R--SQR(  X~2+Y~?+Z“2) 

1420  LET  T2-(  18o*ATN( Y/X) )/3.  14159 

1^*29  IF  INT(I/N1)-(T/M1)sO  THEM  1432 

14  3c  GO  TO  1440 

1431  LET  T2-T2+1S0 

14  3?  PRINT#4:X;n  , " ;Y 

144c  IF  SQR ( X*X+Y* Y ) > 1 THEN  1460 

1441  LET  T4--T0 

144?  LET  Z4 sZ 

1443  LET  WUstf 

1444  LET  K*SQR(X»X+Y*Y) 

1445  IF  T 4 > T 5 + . 0 3 THEM  1452 

1446  IF  K<L  THEN  1456 

1447  LET  D2-D2+ 1 

1443  IF  D2>s2  THEM  1460 

1449  PRINT  2 : T‘3 ; " ";i;"  ";Z5;" 

1450  LET  T5--T5  + 1.4 

1451  GO  TO  146o 

1452  LET  L- 1 

1453  LET  D2-0 

1454  LET  T5-T4 

1455  GO  TO  1460 

1456  LET  LsK 

1457  LET  T5-T4 
1453  LET  Z 5 = 7. 4 

1459  LET  W5sW4 

1460  NEXT  1 

1470 PRINT/* 2:  "FINAL  TIME  ";T0 
143o  PR  I NT/*  3:  TO 
1490  PRINT  // 3 : X , 1 1 

1500PRINT'*?:  "FINAL  Xr  ";X;"VXr  " ; U 
1510  pRIMT *'3 : Y , V 

1520  PRINT/'?:  "FINAL  Ye  " ; y ; " V V e " j V 
1530  PRINT  /'?:  "FINAL  Za";Z;"  VZs";W 
1 540  PRINT  .*>3:7, 

1560  PR  I NT 
1560 END 
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Appendix  IV 


>aputer  Program  - C0NV2 


The  prograa  C0NV2  was  used  to  convert  the  positions  and  velocities 
of  Pluto  in  each  of  the  three  coordinates  that  were  obtained  in  the 
Neptune  frame  (See  pragram  NEPTUNE  in  Appendix  II ) into  positions  and 
velocities  in  the  Sun  center  solar  systea.  These  initial  conditions 
were  given  in  both  the  rotating  and  non-rotating  coordinate  systems. 
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C0NV2 
5 A - 0 

10  PRINT  "INPUT  X,VX,Y,VY,T" 

20  INPUT  X , V 1 , Y , V? , T 
30  LET  0-T*(6.6733629E-4) 

33  LET  Y-Y*3.755E8/4.4999F.1? 

40  LET  XU(X*3.755E8/4.4999E1?+1  ) *C0S(0)-Y*STN(0) 
30  LET  Y1-(X*3«756E3/4.4999E12+1)*SIN(0)+Y*C0S(0) 
60  LET  VUV1*  .736435732-  5. 2331  35303*SIN(0) 

70  LET  V2  = V2*. 736455782+6. 233135303*C0S(0) 

71  LET  V3=V1*C0S(0)-V2*STN(0) 

72  LET  V4=V1*SIN(0)+V2*C0S(0) 

73  LET  X-X1 
73  LET  Y-Y1 

76  LET  V 1 = V 3 

77  V2= V 4 

9o  Bs. 299673033 
100  Y-Y*C0S(B) 

110  V2-V2«COS(B) 

120  Z-Y*SIN(B) 

130  V3=V2»S1N(R) 

191  PRINT  "NON-ROTATING" 

192  PRINT  " X - " ; X ; " V X = " ; V 1 

193  PRINT  " Y +• " ; Y ; " V Y - " ; V 2 

194  PRINT  "Z=" ;Z;"VZ=" ,V3 
193  LET  A = A+1 

196  IF  A- 2 THEN  9^9 
200  V1rV1+6.28318530R*Y 
210  V2=V 2-6. 23318530 R*X 

215  PRINT  "ROTATING" 

216  GO  TO  192 
999  END 
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